Dipolar order by disorder in the classical Heisenberg antiferromagnet on the kagome lattice 



(N 

O 

(N 



I 

o 
o 



o 



Gia-Wei Chern^'^ '^ and R. Moessner^ 
^ Max-Planck-Institut fUr Physik komplexer Systeme, Dresden, Germany 
^ Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico, USA and 
^ Department of Physics, University of Wisconsin, Madison, Wisconsin, USA 
(Dated:'july 20, 2012) 

Ever since the experiments which founded the field of highly frustrated magnetism, the kagome Heisenberg 
antiferromagnet has been the archetypical setting for the study of fluctuation induced exotic ordering. To this day 
the nature of its classical low-temperature state has remained a mystery: the non-linear nature of the fluctuations 
around the exponentially numerous harmonically degenerate ground states has not permitted a controlled theory, 
while its complex energy landscape has precluded numerical simulations at low temperature. Here we present 
an efficient Monte Carlo algorithm which removes the latter obstacle. Our simulations detect a low-temperature 
regime in which correlations saturate at a remarkably small value. Feeding these results into an effective model 
and analyzing the results in the framework of an appropriate field theory implies the presence of long-range 
dipolar spin order with a tripled unit cell. 



The first experiments on the 'kagome bilayer' SCGO 11112] 
triggered a wave of interest in kagome antiferromagnets in 
particular, and frustrated systems in general. A cluster of early 
seminal theoretical papers established kagome magnets 
as model systems for novel ordering phenomena, discussing 
in particular spin liquidity, partial order, disorder-free glassi- 
ness and order by disorder. The excitement persists, not least 
in the quantum realm 10], where there has been much recent 
progress in understanding the ground state for 5* = 1/2 i^. 

Remarkably, for the classical kagome Heisenberg magnet, 
the nature of low-temperature phase has not been established, 
despite the deceptive simplicity of its Hamiltonian, encoding 
only nearest-neighbour interactions of strength J > be- 
tween classical unit-length spins S^: H = Jj2{ij) ' ^j- 
This happens because classical Heisenberg spins do not lend 
themselves at all to the numerical methods applied to 5* = 1/2 
(e.g. jgUTlll). while classical Monte Carlo (e.g. |H|l[l2]) has 
not been able to sample the different local free energy min- 
ima separated by entropic barriers - for the best effort yet, see 
Ref. jlSll . At the same time, analytical approaches have man- 




aged to develop different possible scenarios, described below, 
without being able to choose between them lIsi fTill . 

Here, we present a Monte Carlo algorithm which enables us 
to simulate systems containing over 2000 spins down to a low- 
temperature regime, where we find that correlations no longer 
change as T is lowered further We use the results to deter- 
mine the best parameters in effective models, a stiffness in a 
height model, an effective field theory; as well as a further- 
neighbour coupling in an extended Potts model. In particu- 
lar the latter can then be simulated for over 10^ spins, which 
enables us to verify the critical behavior encapsulated by the 
height model. Thus we identify which ordering scenario ap- 
plies; we obtain spin ordering with a remarkably small or- 
dered moment, which we estimate to be about an order of 
magnitude below fully developed order In the remainder of 
this manuscript, we give an account of our analysis, starting 
with a review of the unusual behaviour of the kagome magnet 
as the temperature is lowered. 

Complex energy landscape: As the kagome magnet is 



FIG. 1. Top: candidate phases of the kagome antiferromagnet (from 
left to right: q — 0; Kosterlitz-Thouless phase; q — \/3 x \/3) 
in the phase diagram of an extended Potts model as defined in the 
text. J2 = corresponds to the pure Potts model on the kagome 
lattice, the unweighted ensemble of coplanar Heisenberg ground 
states. Bottom: cartoon of the effective entropy landscape for ground 
states-soft fluctuations yield minima at co-planar states. Our algo- 
rithm surmounts the entropic free energy barriers, exploring many 
coplanar states. The resulting ensemble average leads to a small but 
nonzero spin order. Inset: transitions between harmonically equiv- 
alent configurations-spins on loops of alternating color are rotated 
such that the angle a spin makes with its local Potts axis is not 
changed as the dark green/blue spins are exchanged for the light ones. 



cooled from high temperatures, it develops short-range order 
like any other magnet. Upon cooling further, its frustrated 
nature asserts itself: at the Curie temperature, J, there is no 
sign of any ordering predicted by mean-field theory, despite 
the fact that the spins in each triangle firmly point at 120° to 
one another. Instead, it is not until much lower T < lO^'^ J 
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that all spins adopt a common plane 115 1 . 

This coplanar ordering is described mathematically by two 
order parameters ifisll : a quadrupolar (also known as nematic 
or coplanar) one for the direction normal to the plane; and an 
octupolar one as follows. If one denotes the angle the spins 
make in the plane with respect to a reference spin by 6, it 
is exp(3i6') which orders. Crucially, the existence of dipolar 
(spin) order, in exp(i0), remains an open question. 

As the ground state is continuously degenerate, the en- 
ergy landscape has many zero-energy directions (flat valleys). 
However, at non-zero temperature, in the free energy land- 
scape - which incorporates entropic effects - the flat valley 
are replaced by robust local minima at each coplanar state. 
The coplanar states are discrete, and can be mapped onto 
the ground states of the antiferromagnetic Potts model on 
the kagome lattice (which is in turn equivalent to the three- 
colouring model of the honeycomb lattice): each spin points 
in one of three directions at mutual angles of 120° with one 
another Labelling these with the colours red, green and blue, 
the coplanar states are those where each triangle hosts a spin 
of each colour Below, we will use the fact that these discrete 
models are exactly soluble lll6ll . e.g. it is known that there are 
exponentially many coplanar states, Nc ^ 1.13^, where N is 
the total number of spins. 

Unusually, the origin of the coplanar ordering is entropic 
- coplanar states have particularly soft fluctuations around 
them: linear spin-wave analysis finds an identical spectrum 
for all of them, including an entire band of zero-energy exci- 
tations js}]. It is thus only when anharmonic fluctuations are 
taken into account that the different coplanar states acquire 
different free energies. No algorithm has so far managed to 
probe these free energy differences, due to the impossibility 
of surmounting the free entropic barriers between coplanar 
ground states sufficiently efficiently. Here we present an al- 
gorithm which achieves this as follows. 

Numerical algorithm: At low temperature, fluctuations 
around an ideal coplanar (Potts) state are small. We can thus 
identify a Potts configuration nearest to the state of the system. 
We then rotate spins on an alternating loop of spins of two 
colours so that these two colours are interchanged. Crucially, 
we have found a way of mapping configurations into each 
other which are equivalent as far as the linearised equations 
of motion are concerned - whether or not to accept a move 
thus depends only on precisely the anharmonic terms respon- 
sible for the low-temperature state selection. Our algorithm 
implements the semiclassical dynamical symmetry identified 
inRef.O 

Firstly, we construct the Potts configuration by assigning 
color R to all spins with a positive projection onto a randomly 
chosen initial spin. Then, the other spins along the resulting 
loops are alternatingly colored G-B. Next, for each loop, we 
evaluate the average orientations of spins colored a =R, 
G or B. Finally, we attempt a move consisting of rotating spins 
G (B) by the angle < (T^, Tq) in opposite directions in the 
nematic plane of the loop; this move is accepted by a standard 
Metropolis condition. To see how the dynamical symmetry is 



implemented, write a spin Sq.; = + i, so that the total 
spin of a triangle i, = Y,a = Eq "^a + La At 
the lowest temperatures, = (the ground state condi- 

tion), and the rotation operation TZrT r = T/j, TZb'^ b = Tf g 
and UgTg = Tb as <(Tb,Tg) = ^. Now, TZrSr,, + 
Ti-BSB.i + T^G^G,i = (^f j^fi ~^f)- The linear equations of 
motion are invariant under this change UtIi . We can thus fully 
equilibrate systems with, depending on boundary conditions, 
up to 2025 spins at temperatures down to T = 10^^ J. Lower 
temperatures are achievable but, as we discuss next, this is not 
necessary. 




FIG. 2. Spin correlation functions C3(r) = (S(r) • S(0)) along the 
nearest-neighbor directions at various temperatures. Also shown for 
comparison is the case of Potts model. The simulations were done 
on a lattice with ~ 9 x 15^ spins with open boundary conditions. 
Inset: the correlations acquire a stable low-temperature value below 
T = 10"^ J. 

Correlations saturate at low temperature: There is a lim- 
iting low-temperature regime where the correlations cease to 
change as T is lowered further below T ^ 10^'^ J. This is 
shown in Fig. |2]for spin-correlations in real space for a system 
of = 2025 spins with open boundary conditions. This im- 
plies that the effective weights differentiating between differ- 
ent coplanar states become temperature independent, as fol- 
lows. The fluctuations selecting an ordered state are subdom- 
inant compared to those which select the nematic order: the 
modes which cost no energy are dressed by the finite-energy 
ones when anharmonic effects are taken into account. To low- 
est order, cubic (quartic) terms appear to second (first) order 
in the corresponding perturbation theory. The corresponding 
contributions are suppressed by a factor of T compared to the 
leading ones |14:]. Temperature therefore drops out of the re- 
sulting effective weights. The low-temperature limit is there- 
fore smooth. Our numerics establishes that this holds for the 
full problem, and not just to leading order in a perturbation 
theory which may or may not converge on account of the large 
number of zero modes. We will use this property extensively 
below. 

Enhanced correlations at low T: In Fig. [3] we present the 
correlators necessary for the demonstration of order in this 
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regime of low temperatures. We focus on the order param- 
eter corresponding to a tripling of the unit cell, m^, cor- 
responding to the right state depicted in the phase diagram 
Fig-IH Shown for comparison are the correlations of the (un- 
weighted) average over the coplanar Potts states, which do 
not take into account the anharmonic fluctuations responsi- 
ble for the low-temperature regime under consideration. The 
771^3 deviates upwards from the decay of the Potts correla- 
tions, while the decay of mo, corresponding to the left state in 
Fig-IU is even faster than at higher temperatures (not shown). 
Next, we analyse this data in detail with the aid of effective 
field theories and an effective Hamiltonian. 
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FIG. 3. Top: spin correlations of the Heisenberg model atT — 10~ 
and the extended Potts model with J2 = 0.019 on a lattice with 396 
spins. Bottom: order parameter as function of system size for 
different values of J2. The solid red points for small systems are 
from the Heisenberg simulations. Inset: scaling collapse of data for 
L > 60. 

Height model analysis: Historically, the question about the 
nature of the ordering has sometimes been phrased as an al- 
ternative between ordering in either q = or in a 
structure as these are frequent outcomes of perturbations to 
the classical Heisenberg model. However, the actual distinc- 
tion has to be made between a power-law correlated state and 
an ordered state far away from the q = state, as suggested 
by the phase diagram in Fig. [T] This was noted long ago by 
Huse and Rutenberg |3|], based on Baxter's observation that 
the Potts model, i.e. the average over coplanar states without 
entropic weighting, is at a critical point in the limit of T 
1 1611 . To settle the question whether the kagome magnet is 



in turn ordered, one therefore needs to establish whether the 
appropriate tuning parameter - a stiffness - increases rather 
than decreasing or remaining the same: for decreasing stiff- 
ness, the system would move away from the critical point into 
the KT phase, whereas for an increase, ordering ensues. 

Given our algorithm can access the low-temperature corre- 
lations, we readily find that the stiffness has increased. This 
is evident already in the real-space correlations, which are 
significantly enhanced for Heisenberg over Potts (Fig. |2]i. 
More directly the stiffness can be extracted by considering 
long-wavelength fluctuations of a so-called height field h(r) 
defined at each hexagon of the lattice. The height differ- 
ence between two neighboring hexagons (labelled by /, J) is 
h/ — hj = Yf J, where IJ denotes the site shared by the 
two hexagons jlSll . The ratio of the stiffnesses is given by 
limg^o©!?^) = |hpotts(g)|V|hHcisonborg('7)P, as shown in 
Fig. |4] This quantity, clearly bigger than one, does not even 
appear to reach a limiting value, indicating not only enhanced 
correlations but even the onset of ordering. The nature of the 
ordering is indicated by the value of the average height (which 
is a periodic variable and hence best plotted over its hexagonal 
height 'Brioullin zone') lIlSll . The maximum (minimum) at its 
K{T) point correspond to the q = \/3 x V3{q = 0) states, re- 
spectively, and the difference between these extrema is greatly 
enhanced over the unweighted Potts model (not shown). 

This establishes that the Heisenberg antiferromagnet is 
more ordered than the Potts model which describes the inter- 
mediate, coplanar, temperature regime: the magnet finds itself 
on the ordered side of the transition. 

Effective Hamiltonian: The weakness of the resulting or- 
der turns out to be a blessing in disguise: it allows us to use 
an effective low-temperature Hamiltonian for the Heisenberg 
magnet which is only a small perturbation of the antiferro- 
magnetic Potts model. This is demonstrated in Fig. [3] where 
we find quantitative agreement for the spin correlations at all 
distances between Heisenberg and an extended Potts model 
obtained by merely adding a weak, J2 = 0.019, second- 
neighbour ferromagnetic coupling. This interaction is nor- 
malized such that (un)equal Potts spins (lose) gain an 'en- 
tropy' J2, i.e. the concomitant Boltzmann factors are of the 
temperature-independent form cxp(±, 72), as explained above. 
Of course, there may be small additional terms which we can- 
not accurately fit to the data of the present simulations but the 
agreement between the correlations from the two models is 
remarkable given the simplicity of the fit, and we believe this 
is a basis for a good semi-quantitative approximation scheme. 

With this in hand, we can now simulate the effective Potts 
model to get a better handle on very large system sizes. Fig. [3] 
displays the size of the for different values of J2 for 
systems of up to over N = 10® sites. While does not 
saturate for even those large systems, it never approaches a 
limiting algebraic behavior either and instead always bends 
upwards. Even though ordering is not directly detectable in 
an unambiguous fashion, we can use information from the ex- 
act solution, which leads us to expect a transition at J2 = 0. 
Assuming a scaling form, appropriate for a Kosterlitz- 
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induced order very much smaller than the maximal possible 
moment, which is rather unusual for classical systems. These 
tend to be characterized by robust order, such as the saturated 
ground state order of the frustrated triangular lattice Heisen- 
berg antiferromagnet. Like in the case of the kagome S = 1/2 
quantum magnet, our classical phase is also extremely fragile, 
destabilized already by additional interactions of strength of 
only a few percent of the leading one. 

Note added: As we were concluding this manuscript. 



Ref. |23| appeared, which also considers the low-temperature 
behavior of the classical kagome Heisenberg magnet. Their 
main interest is the dynamics but they do also present an algo- 
rithm accessing the low-temperature regime. 
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FIG. 4. Top: A ratio Q{q) > 1 for reflects the in- 

creased stiffness of the Heisenberg model compared to the Potts 
model, and the resultant ordering at low temperature. Bottom: Rela- 
tive frequency of the averaged height vector for Heisenberg spins at 
T = lO""* J for a lattice with 729 spins and open boundary condi- 
tions. 



Thouless transition lfl9ll . 
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we obtain data collapse for m^. Fig. [5] a, the one fitting 
parameter left indicates a small but non-zero ordered moment 
about an order of magnitude below the full moment. 

Conclusions and outlook: We have analyzed the asymptotic 
low-temperature regime of the classical kagome Heisenberg 
antiferromagnet with a combination of analytical and numeri- 
cal approaches. We find dipolar spin order with a tripled unit 
cell and a very small ordered moment. It will be interesting to 
see how this ordering disappears as the temperature is raised. 
More broadly, it is desirable to study the effective model in 
considerably more detail as a function of J2, in particular with 
view to the question of whether there are any additional tran- 
sitions, e.g. involving the chiralitv 12011 . the other marginal op- 
erator of the critical Potts model 11611 . Additional interesting 
models to examine would be the ordering of analogous mag- 
nets on the pyrochlore lattice |f2ltl, for some work on this see 
Ref. II22II . or the fate of the hyperkagome magnet in the limit 
of the lowest temperatures 11311 . 

Overall, this kagome magnet provides a most striking 
instance of a classical system with nonlinear-fluctuation- 
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